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Abstract 

We study relations between the ground-state energy of a quantum graph Hamiltonian with attractive S coupling at the 
vertices and the graph geometry. We derive a necessary and sufficient condition under which the energy increases 
with the increase of graph edge lengths. We show that this is always the case if the graph has no branchings while 
both energy increase and decrease are possible for graphs with a more complicated topology. 
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1. Introduction 

Quantum graphs proved themselves to be a class of systems offering numerous problems interesting from both the 
physical and mathematical point of view; we refer to the proceedings volume [Q]] for an extensive bibliography. In this 
Letter we address the question about relations between the ground-state energy of such a Hamiltonian and geometric 
properties of the underlying graph, in particular, the lengths of its edges. 

A motivation to study this kind of problem is twofold. On the physics side it is, of course, the importance of 
the ground state as the one to which the system tends to relax when it loses energy due to an interaction with the 
environment. Since quantum graphs model various real physical systems it is natural to ask about the geometric 
configurations which are energetically the most favourable. At the same time, mathematically the problem represents 
a natural extension of the usual spectral-geometry studies of the relations between spectral properties of differential 
operators and geometry of the manifolds supporting them. 

We restrict here our attention to graphs with a finite number of edges, some of which may be semi-infinite, and 
an attractive 6 coupling at the vertices, assuming that the motion at the graph edges away from the vertices is free. 
Such systems have always a nontrivial negative spectrum with a well-defined ground state; we will ask how the 
corresponding eigenvalue depends on the finite-edge lengths. First we analyze the case of n attractive 5 interactions 
on the line which can be regarded as a simple chain graph. We will prove that the ground-state energy moves up with 
increasing distances between the 8 potentials in two different ways, by means of a Neumann bracketing and by using 
the well-known explicit form of such a Hamiltonian resolvent. 

After that we will pass to general quantum graphs of the described class. We will show that in such a case the 
dependence on the edge length is more complicated and its sign is uniquely determined by the form of the ground-state 
eigenfunction on the particular edge. As long as the graph is a chain we have the monotonicity described above. On 
the other hand, we will give an example showing that once the graph has at least one nontrivial branching, i.e. a vertex 
of degree exceeding two, it is possible that the ground-state energy decreases with the increasing edge lengths. 

Before proceeding further let us note that relations between quantum graph eigenvalues and edge lengths have 
been discussed also in other contexts. In particular, Friedlander [2J] derived a lower bound on higher eigenvalues for 
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finite graphs in terms of the total graph size. On the other hand, Berkolaiko and Kuchment j3s] studied general relations 
between the point spectrum and the set of edge lengths and coupling constants. 

2. A warm-up: 6 interactions on a line 

Consider first a particle on a line with a finite number of (5-interactions the Hamiltonian of which can be formally 
written as — jjjp + Y!)=\ oijS(x - yj). Following |4j] we denote this operator as -A a j where a := {ai, . . . , a„) and 
Y :- {yi, . . . ,y n ). We suppose that all the points yj are mutually distinct and the interactions are attractive, aj < 0, 
j = l,...,n. Under this assumption the continuous spectrum of -A ff y covers the positive halfline and the discrete 
spectrum in the negative part of the axis is non-empty, in particular, there is a ground-state eigenvalue Aq < with a 
strictly positive eigenfunction (/^fl; we ask how does Aq depend on the geometry of the set Y. 

One can conjecture that the ground-state energy decreases when the point interactions are closer to each other. 
First we prove this claim under an additional assumption. 

Proposition 2.1. Consider sets Y\, Y% of the same cardinality such that yj t i < yj t 2 < . . . < yp, , j — 1, 2. Let there 
be an i such that y2j — yij for I — 1, . . . , i and y2j — yij + t] for I = i + 1, . . .,n. Suppose further the ground-state 
eigenfunction of the — A a j i satisfies (^Cyi,;+) < and tft' Q (yij+i— ) > 0. Then we have min<x(— A^y,) < mincr^-A^y.,) 
for any tj > 0. 

Proof: Since i^o is positive and satisfies t/^' = -Aq^o between the point interaction sites, the function is convex; 
by the assumption there is then a point xq e (yi.^Vi.j+i) such that 4t' (xo) = 0. Consider now the operator -A a y, 
which acts as -A tr ,y, with the additional splitting^ Neumann condition at the point xo; it is obvious that the two 
operators have the same ground state. Such a Neumann condition separates the two halflines, hence -A ff y, can be 
written as -A^ Yi © —A r a Fj . Consider now the operator -A tr ,y, := -A' v Yi ffi -A^ ffi -A' a Yi where the added operator 
is the Neumann Laplacian on L 2 (0, rf)\ it is clear that the latter does not contribute to the negative spectrum, hence 
mino^-A^y,) = mincK-A^y,). Furthermore, -A a y, is obviously unitarily equivalent to -A ff y, with added splitting 
Neumann conditions at the points x = jco, xq + 77, hence the sought result follows from Neumann bracketing ||5l 
Sec. XIII.15]. □ 

It is not difficult to see that the assumption about the derivative signs is satisfied if -a-,, -a,+i are large enough or, 
which is the same by scaling, the distance yi+\ - y, is large enough. However, we can make a stronger claim without 
imposing restrictions on the ground-state eigenfunction derivatives. 

Theorem 2.2. Suppose again that #Y\ — #Y2 and aj < for all j. Let further yij — y\j < y^i — y2j hold for all i, j 
and yij — yij < y2j — y2jfor at least one pair of i, j, then we have min <x(— A^y, ) < min cr(— A a y 2 ). 

Proof: We employ Krein's formula |0, Sec. II. 2. 1] which makes it possible to reduce the spectral problem at energy k 2 
to solution of the secular equation, detF ffi y(fc) = 0, where 

[T a , Y (k)]jf = + Gkiyj - X/')]J/=i 

and Gkiyj - y/) = jre lk \ yi ~ y >'\ is the free resolvent kernel. Writing conventionally k = Ik with k > 0, we have to 
investigate the lowest eigenvalue of V a j(iK) which is, of course, given by 

p.o(a, Y; k) = min (c, F a y(Ik)c) 

\c\ = l 

with the minimum taken over all c e C" with |c| = 1 . It is easy to see that po(a, Y; k) > for all k large enough; 
the ground state energy -k 2 corresponds to the highest value of k such that fio(a, Y\ k) = 0. Since [T a j(iK)]ij = 



'See d Thm. II.2.1.3], and also TheoremHHbelow. 
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condition is sometimes used as a synonym for Kirchhoff coupling in quantum graphs, we say "splitting" to stress that the functions from the domain 
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—5ija i 1 - y e where ^ = |j; - the quantity to be minimized is explicitly 

(c r a , Y (i K )c) = J] |c,-| 2 1-- - -J - 2 2 J] Re 

Next we notice that the eigenfunction corresponding to the ground state, i.e. c for which the minimum is reached can 
be chosen strictly positive; we write symbolically c > meaning c, > 0, i = 1, . . . ,n. This follows from the fact 
that the semigroup {er tr "- T(fK ' : t > 0} is positivity improving, as a consequence of strict negativity of the off-diagonal 
elements of T a j(iK) — cf. fl, Sec. XIII. 12 and Problem XIII. 97. This means, in particular, that we have 

uo(a, Y; k) = min (c, F a y(//c)c) 

|c|=l,c>0 

Take now two configurations, (a, Y) and (a, Y) such that £y < Ijj and the inequality is strict for at least one pair (7, j). 
For any fixed c > we then have (c, r Q . j y(i/e)e) < (c, F a y(iK)cj, and consequently, taking a minimum overs all such c's 
we get 

/i (a, Y; k) < p. (a, Y; k) 

for any k > with the obvious consequence for the ground state of -A a j', the sharp inequality in the last formula 
holds due to the fact that there is a c for which the minimum is attained. □ 

Remark 2.3. The argument used above can be extended to other situation. Take for instance, point interactions on a 
loop, in other words, on a finite interval with periodic boundary conditions. The corresponding Green's function is 

G iK (x, y) = — — — , \x-y\<-€, 

2k sinh kc 2 

where 6 is the length of the loop. Writing the corresponding secular equation we find that expanding the loop without 
reducing the distances between the neighbouring point interaction sites means moving the ground-state energy up. 

3. Quantum graphs: setting the problem 

After this preliminary let us pass to a more general situation when the particle lives on a graph and the attractive 
point interaction represent couplings at the graph vertices. Consider a graph T consisting of a set of vertices *V = 
{Xj : j e /}, a set of finite edges X = {-Cj n ■ (j, n) e l£ <z I x 1} where £j n is the the edged connecting the vertices Xj 
and X„, and a set of infinite edges Xco = {-Geo : k e Iq) attached to them. We regard it as a configuration space of a 
quantum system with the Hilbert space 

<H = L 2 ([0, /,■]) e L 2 ([0, oo)) . 

jel z kel c 

the elements of which can be written as columns i/r = {(tffj„ : -Cj„ e -C), {if/kco ■ -Ckoo e -Coo}) T ■ We consider the 
dynamics governed by a Hamiltonian which acts as -d 2 /dx 2 on each edge. In order to make it a self-adjoint operator, 
in general boundary conditions 

(Uj-iy¥j + i(Uj+I)Tj = (3.1) 

with unitary matrices Uj have to be imposed at the vertices Xj, where TJ; and *F. are vectors of the functional values 
and of the (outward) derivatives at the particular vertex, respectively la, LZl |8[]. In other words, the domain of the 
Hamiltonian consists of all functions in W 2,2 (£ ffi -Coo) which satisfy the conditions ( 13. U . We will be interested in the 
following particular class: 

• the internal part of the graphs is finite and so is the number of external edges, #/_£ < oo and #Iq < oo 



3 Without loss of generality we may suppose that each pair of vertices is connected by a single edge; in the opposite case we add extra vertices 
of degree two to the "superfluous" edges and impose Kirchhoff conditions there. 
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• the coupling at each vertex is of 6 type in terminology of |9|], i.e. Uj = —f^-^T - 1, where nj is the degree of the 
vertex Xj and J' is the matrix having all the entries equal to one. Explicitly the coupling conditions ( 13.11 ) then 
become 

<A/,(0) = -MO) =: ^(0) , j e I, k = 1, . . . , nj , J] W> = a ^' (0) ' (3 ' 2) 

i=l 

where each edge emanating from Xj is parametrized in such way that x — corresponds to the vertex 

• for "free endpoints", or vertices of degree one, parametrized by xj = I j, this in particular means the Robin 
condition, iff'j(lj) + aji//j(lj) = 

• all the couplings involved are non-repulsive, ay < for all j e /, and at least one of them is attractive, aj < 
for some jo e / 

In such a case it is not difficult to express the quadratic form associated with the quantum-graph Hamiltonian H: 
it is given b>0 

?m = 2 fV«l 2 d* + X f l^(*)| 2 dx + £ ai |^(0)| 2 , (3.3) 

jel £ J° kel c iel 

where ifrj, fa are components of the wave function *P on the internal and external edges, respectively, and ^,-(0) are 
the values at the vertices. The domain of the form consists of Lr functions which are W 1,2 on the graph edges and 
continuous at the vertices. 

Proposition 3.1. inf cr(H) < holds under the stated assumptions. 

Proof: If Iq — we take a constant function, VP = c on T which belongs to the form domain because T has then a finite 
length; we get q\^¥\ < a ;o |c| 2 . On the other hand, if Iq + 0, we take *P equal to c on the internal part of the graph and 
to iffk(x) = c e~ KX on each external semi-infinite edge. The integrals over the internal edges vanish as before and those 
over external ones are easily evaluated; we get 

q{V] < \a Jo + i*#/ c J |c| 2 
which can be made negative by choosing k small enough. □ 

Theorem 3.2. In addition, let Y be connected, then the bottom of the spectrum Aq — inf cr(H) is a simple isolated 
eigenvalue. The corresponding eigenfunction T' ' can be chosen strictly positive on V being convex on each edge. 

Proof: Consider a disjoint graph with all the vertex couplings changed to Dirichlet conditions. In such a case the 
spectrum is positive; it is discrete if Iq — and equal to R+ otherwise. By Krein's formula (10, Proposition 2.3], the 
original operator differs from the Dirichlet decoupled one by a finite-rank perturbation in the resolvent, hence their 
essential spectra are the same by Weyl's theorem and the negative spectrum of H may consist at most of a finite number 
of eigenvalues of finite multiplicity; by the previous proposition it is nonempty and the ground-state eigenvalue exists. 



The ground state positivity follows, e.g., from a quantum-graph modification of the Courant theorem 111 111 . The 
eigenfunction being positive and its component i/A ' at the y'th edge twice differentiable away of the vertices, we have 

(<// 0) )" = -Aoil/j > 0, which means the convexity. □ 

In fact, one can say more about the ground-state eigenfunction because the corresponding Schrodinger equation 
can be solved explicitly. Writing the spectral threshold as Aq = -k 2 we see that the eigenfunction component on each 



A meticulous reader might notice that numberings of the functions on the edges differ; sometimes it is practical to number the edges, sometimes 
vertices at their endpoints, or edges sprouting from a given one. We are sure that this can cause no misunderstandings. 
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edge is a linear combination of e n and e KX . Since we are free to choose the edge orientation, each component has 
one of the following three forms, 

Icj cosh k(x + dj) , dj e R 

cje****^, dj&R (3.4) 

Cj sinh k(x + c/y) , x + dj > 

where cy is a positive constant. For further purposes we introduce edge index 

+ 1 ... ^\x) = Cj cosh k(x + dj) 

... ifrJ\x) = c J e tK( - x * d i ) (3.5) 
-1 



<x ; := < 



tpj\x) = Cj sinh k(x + dj) 



4. Monotonicity proof by a scaling argument 



From now on we consider connected graphs only which we can do without loss of generality, since otherwise we 
deal with each connected component separately. By Theorem l3.2l the graph's Hamiltonian H then has a simple ground 
state with positive eigenfunction. Using the above definition, we can compare graphs with the same index structure. 

Given F and f with the same topology differing possibly by inner edge lengths, we consider the family of inter- 
polating graphs having the length of the j'th edge in the closed interval with lj and lj as endpointfl We say that the 
graphs F and f belongs to the same ground-state class if the edge indices of the graphs edges remain the same for this 
whole family. Equipped with this notion we can make the following claim. 

Theorem 4.1. Under the stated assumptions, consider graphs F and F of the same ground-state class. Let H and H be 
the corresponding Hamiltonians with the same couplings in the respective vertices, and Aq and Ao the corresponding 
ground-state eigenvalues. Suppose that o-jlj < o-jlj holds all j 6 /r such that \o-j\ — 1 and lj — lj if o~j — 0, then 
Aq < Aq; the inequality is sharp if o-jlj < o-jlj holds for at least one j € /_£. 

Proof: It is obviously sufficient to compare graphs differing just by the length of a single inner edge corresponding to 
a fixed index value j e Iq with \o~j\ = 1, and it is enough to prove the claim locally. We choose a finite-length segment 
J = [a,b] in the interior of the y'th edge and write F as the union of J and Fj := F \ J. Without loss of generality 
we may choose J in such a way that b - a > lj - lj if <x,- = 1 and b - a < I j if o~j = -1. Then f can be written as 
FjU J where J is obtained by scaling of J with the factor £ := being less than one in case of a shrinking edge 

and larger than one otherwise. In order to prove the desired result we have to find a function *F e L 2 (F) such that the 
Rayleigh quotient on the tilded graph satisfies 

(4.1) 

imp v ; 

We construct such a trial function *F (0) in the following way: we put ^ (0 \x) = ^^{x) for x e Fj while the y'th 
component on / is obtained by scaling 

trf\a + fy)=i{,f(a + y) for 0<y<|/|; (4.2) 

in order to prove ( 14.11 ) we have to choose £ < 1 iff o~j = 1 and vice versa. The Rayleigh quotient for the function ( 14.2b 
can be easily rewritten in a natural notation as 

= — =: /(f) . (4.3) 



5 This family corresponds to a closed parallelopiped in the natural parameter space (0, oo)"'x of our problem given by the interior edge lengths. 
If at least one of the edges has the same length in both the T and f" the said parallelopiped is degenerate and we regard it instead as a subset in the 
reduced parameter space referring to the changing edge lenghts. The eigenvalues are analytic functions on the interior of such a set - cf. 0]. 
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where 

a-qrjW®], b := Juf )'(x)\ 2 dx, 

and c, d are the parts of the squared norm of ^F' ' corresponding to Fj and /, respectively. It is enough to check that 
<Tjf'{\) = -<Tj(bc + 2bd + ad)(c + d)~ 2 > 0. Choosing the ground-state eigenfunction *P (0) conventionally with the 
norm equal to one, we have c + d = 1 and a + b — Aq, hence the property to be checked is -o-j(Aod + b) > 0, or more 
explicitly 

-^(^oll< ) H 2 J + ll< ) )'ll 2 7 )>0. 

Using Aq = -k 2 we find for o~j = 1 

J \{il/f)'(x)\ 2 Ax = c 2 k 2 J^( sinh KX f Ax < c)* 2 J^(cosh kx) 2 dx 
= -A JyP(x)\ 2 dx, 

and the opposite inequality for o~j = -1 where the roles of hyperbolic sine and cosine are interchanged. Hence 
inequality (14. j} is satisfied for for £ < 1 if crj = 1 and £ > 1 if o~j = — 1, provided |1 — £| is small enough, which is 
what we have set out to prove. □ 

Remark 4.2. The case o~j = is nontrivial; the critical case in the example given in the next section shows that such 
a possibility is not excluded in case of finite edges and the ground-state energy can be independent of length changes 
of such edges, however, a more subtle analysis is needed to treat such situations in general. 



5. Discussion 

Let us now ask what consequences can one derive from our main result given in Theorem 14. II First we notice 
that the graphs without branchings belong all to the same class and one is able to extend to them conclusions of 
Theorem l2.2l and Remark l2~3l Specifically, we can make the following claim. 

Corrolary 5.1. In the setting ofTheorem \4.1\ suppose that the graph T has no branchings, i.e. the degree of no vertex 
exceeds two. Then the index of any edge is non-negative being equal to one for any internal edge. Consequently, a 
length increase of any internal edge moves the ground-state energy up. 

Proof: By assumption a graph without branchings is a chain of edges, either closed into a loop or open; in view of 
Remark 12.31 we can consider only the latter possibility. It is obvious that it is not possible that all the edges have 
negative index. If the first or the last one are semi-infinite, their index must be zero; if all the edges are finite the 
attractive (5-coupling would require that sgn if/ remains unchanged over the whole chain and the non-repulsive Robin 
condition on one of the endpoints could not be satisfied. The question is whether one can have a sinh-type solution 
at some position within the chain. In such a case there would be a vertex in which wavefunction components with 
different indices have to match. Let us parametrize the chain by a single variable x choosing x = for the vertex 
in question. Suppose that the (non-normalized) ground-state eigenfunction equals if/j{x) = cosh/f(c/i - x) for x < 
and \j/j + \{x) = c sinh^^ + x) for x > 0. By assumption they are coupled by an attractive 8 interaction, hence c 
is determined by the continuity requirement and ^ +1 (0+) - t^(O-) must be negative; recall that the ground-state 
eigenfunction is positive. However, this expression equals +k cosh K{d\ + do) I sinh/c^, hence the needed match is 
impossible for a sinh solution decreasing towards the vertex. The same is true for the opposite order of the two 
solutions, and similarly one can check that a negative-index edge cannot neighbour with a semi-infinite one. □ 

On the other hand, for graphs with a more complicated topology the analogous claim is no longer true. We will 
illustrate it on a simple example of a star graph with mirror symmetry sketched on Fig. Q] We have plotted here the 
ground-state energy — in the logarithmic scale to make the effect more visible — as a function of the edge length Lo 
and the coupling constant a in the central vertex. We see two regimes here. For weak attractive coupling, a cr it < a <0 
where or cl -i t » -1.09088, the ground-state energy decreases with increasing L2 while the opposite is true if a < o- CI i t ; at 
the critical value the energy is independent of Lo and the solution on the "axial" edge is a pure exponential. 
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Figure 1 : The ground-state energy of sketched star graph as a function of L2 and a. 

The reason why this happens in the example is obvious. The mirror symmetry allows us to decompose the problem 
into a symmetric part, where the ground state is to be sought, and antisymmetric one which reduced trivially to the 
Dirichlet problem on a single interval. Using the notation from the proof of Corollary 15.11 the left-hand side of the 
derivative condition in the symmetric part equals if/' 2 (0+) - 21/^(0-), hence the argument used there no longer applies. 
On the other hand, it is not difficult to construct examples without a symmetry in which we have different regimes; an 
open question is whether one can find a general regime characterization for an arbitrary branching graph. 

Let us finally recall that our main result, Theorem l4.ll holds for graphs belonging to the same ground-state class. 
We know from yl that the notion is not empty: it follows from eigenfunction dependence on the edge lengths that 
a given graph has a neighbourhood in the parameter space where the indices do not change. The question about 
existence of different classes for graphs of the same topology is open and interesting. 
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